Two-body scattering without angular-momentum decomposition 
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Two-body scattering is studied by solving the Lippmann-Schwinger equation in momentum space 
without angular-momentum decomposition for a local spin dependent short range interaction plus 
Coulomb. The screening and renormalization approach is employed to treat the Coulomb interaction. 
Benchmark calculations are performed by comparing our procedure with partial-wave calculations 



in configuration space for p- 
potential model. 



"Be, p-^^O and ^^C-^°Be elastic scattering, using a simple optical 
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I. INTRODUCTION 

The aim of the present work is to solve the two-body 
Lippmann-Schwinger equation without partial-wave de- 
composition for a local short range interaction plus 
Coulomb. This is a first step towards the ultimate goal of 
solving exact three-body equations without partial-wave 
decomposition as a means to describe complex nuclear 
reactions where three-body degrees of freedom play a sig- 
nificant role. 

The inclusion of the long range Coulomb force be- 
tween charged particles of equal sign has become pos- 
sible, in recent years, through a novel implementation 
of the method of Coulomb screening and renormaliza- 
tion [l], [3| in the framework of Alt, Grassberger and 
Sandhas (AGS) exact three- 01 and four-body in- 
tegral equations, leading to fully converged results for 
three- and four-nucleon scattering [7', and for di- 
rect nuclear reactions dominated by three-body degrees 
of freedom 0, E [U, [ll • 

In all these calculations [1, 0, 0, S, B E, [HI, [H the 
equations were solved using the partial-wave decomposi- 
tion of the multivariable integral equations and taking as 
many partial waves as needed for convergence of the ob- 
servables. Although we get fully converged results, at in- 
termediate energies the partial wave expansion converges 
very slowly and may get unstable if we increase the en- 
ergy beyond the values we used or if we address reactions 
with two heavier nuclei such as ^^C-^^Be where ^^C and 
^•^Be are considered as inert cores. In this case, even 
at moderate energies the relative ^^C-^^Be wave length 
is so small that a very large number of partial waves is 
required for convergence. 

Therefore, in order to get exact three-body results at 
higher energies or involving the collision of two heavy nu- 
clei one may need to develop calculations without partial- 
wave decomposition. The group at Ohio University has 
already made progress in this direction for two-body 
[13. . .14 1 and three-body scattering (isl . [l6| . following 
an earlier work by Belyaev et al. However they do 
not include the exact treatment of Coulomb in their cal- 



culations, though in their early works [H, [l^, they 
introduced Coulomb in an approximate way in the con- 
text of a multiple scattering framework. 

In the present manuscript, we show results for the solu- 
tion of the two-body Lippmann-Schwinger equation for 
p-^^Be, and ^^C-^°Be elastic scattering at inter- 

mediate energies. The novelty of the present work vis- 
a-vis the Ohio group two-body calculations is that we 
include Coulomb between two charged cores of mass Ai 
and atomic number Zi together with standard optical 
potentials that include a central plus a spin-orbit inter- 
action. 

This paper is structured as following. In Sec. II we 
present the formalism of the Lippmann-Schwinger equa- 
tion for central and the spin-orbit potentials. In Sec. Ill 
we summarize the Coulomb treatment. In Sec. IV we 
show the results for different reactions. In Sec. V we 
summarize and give conclusions. Finally in Appendix A 
we present the analytical expressions for the Fourier 
transform of the potentials used in this work and in 
Appendix B we explain the method used to solve the 
Lippmann-Schwinger equation. 



II. LIPPMANN-SCHWINGER EQUATION 

The Lippmann-Schwinger equation for the scattering 
of two particles is 



T^V + VGqT, 



(1) 



where V is the two-body potential between the parti- 
cles. Go = {Z — Ho)~^ the free two-body resolvent, and 
T the transition operator. The matrix element of the 
transition amplitude in momentum space, T(q', q, Z) = 
(q'|T(Z)|q), satisfies the integral equation 

T(q',q,Z) = T/(q',q) (2) 
+ / d^q"V{ci\q") \,^,„ r(q",q,Z). 

J Zl — 7\ 

Here, q is the relative wave vector, // is the reduced mass 
of the two particles and Z the appropriate energy. 
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Central interaction 



First we consider the case of a central interaction only 
where we follow the same procedure as Elster et al. [l3 |. 
In this case, the matrix elements in momentum space of 
the potential and the transition amplitude, V(q',q) and 
T(q', q, Z) respectively, are scalar functions 

V'(q',q) = V{q',q,q' -(i), (3) 
T(q',q,Z) = r(g',g,q'.q,^). (4) 

Therefore Eq. ([2| can be expressed as follows 

T(g', 9, x', Z) = V{q', q, x') + dq" q"^ j ^ dx" 



X I dip"V{q\q",y)--^Tiq",q,x",Z),i5) 



2n 



where a;' = q' • q, x" — q" • q, and y = q" • q'. We take 
the incoming wave vector q in the direction of the z-axis 
and the arbitrary azimuthal angle ip' for q' is chosen to 
be zero. We can then express y through x' and x" as 



y — x'x" + vT— a/^yl 



x"^ cos (f" . 



(6) 



Defining 

viq\q",x',x")^ r d^"V{q',q",yl (7) 
Jo 

the integral Eq. ^ becomes 

T{q\ q, x\ Z) = —v{q', q, x' , 1) + dq"q"' (8) 



X / dx"v{q',q",x',x") 



Z - 



R2, 



—T{q",q,x",Z), 



2ti 



leading to a two-dimensional integral equation in the off- 
shell wave vector q" and the cosine of the scattering angle 
x". 



B. Spin-orbit interaction 

Next, we consider the case in which we have a spin- 
orbit interaction. In this case, Eq. ^ becomes a set of 
coupled equations 

Ts'A',.A(q',q,^) = T4'A',.A(q',q) (9) 

+ E / rf-\"^4'V,."A"(q',q") 

s" X" 

-r."A".sA(q",q,^), 



where s is the spin of the system and A its projection in 
the z-axis. The spin-orbit potential, commonly expressed 
as 



is not central anymore. The a ■ 1 term introduces a de- 
pendence on the azimuthal angle (/?' which makes Eq. ^ 
not valid. Nevertheless, it is possible to reduce Eq. ^ to 
a two-variable integral equation. From the partial-wave 
analysis of the T-matrix 

T,,v,.A(q',q,Z) = E E Y.YL'M'm (11) 

JMj L'M' LM 

X {L'M's'\'\JMj)Ty''''\q', q, Z){LMs\\JAIj)Ylj,j{ci), 

assuming, as before, the initial wave vector q to be 
along the z-axis, it follows that the (^'-dependence of 
Ts'\',s\{(l' ,1, Z) is determined by the spherical harmon- 
ics rL'Af'(q') in terms of e**^'"^' with fixed M' = X - X' 
because M = 0. Therefore the T-matrix can be written 
in factorized form as 

Ts'X'Mq.', q, z) = e'^^-^">'^'%'yMq\ z) (12) 

where Tg'y ^s\{q' , q,x' , z) is the solution of a set of two- 
variable integral equations 



Ts'x'^s\{q',q,x',Z) = t^-Ws'A'.saI?'. 9- 2^'. 1) 
zvr 



(13) 



/ JJIJI2 I j„//„A 



s"A" 



dq"q" 



dx"v^,^,^,„y,{q',q\x\x") 



Z 



-772%"x",s\{q" ,q,x" ,Z) 



2/i 



Here v^,y ^„y,{q' ,q" ,x' ,x") includes the phase from the 
T-matrix and is chosen to be zero 

4a',."A" («',?", 2:', x") (14) 



2tt 



d^" e'(^-^")'^'V,,v,,.A.(q',q")l^' 



Vso{^) = Vso{r) a ■ 1, 



(10) 



¥''=0- 



As discussed in Ref. [H, [15I . [2H these calculations 
are time consuming if the potential is given in config- 
uration space and the transform to momentum space 
(Fourier transform) is performed numerically. Therefore 
in the present work we develop in Appendix [A] analytic 
Fourier transforms for Woods-Saxon interactions, cen- 
tral, surface, and spin-orbit, together with the screened 
Coulomb interaction. In Appendix [B] we outline the nu- 
merical procedure we follow to solve the integral equa- 
tions. 



III. TREATMENT OF COULOMB 
INTERACTION 

The inclusion of the Coulomb interaction in momen- 
tum space is a very complicated task due to its l/g^ 
singularity which together with the G^iZ) singularity 
renders the kernel of the Lippmann-Schwinger equation 
non-compact. This fact has been a handicap for perform- 
ing momentum space scattering calculations involving 
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charged particles. Over the years several methods have 
been proposed to overcome the Coulomb singularity that 
introduce a cutoff parameter: the pioneer work was done 
by Vincent and Phatak [l^l which is, in principle, ex- 
act and works well for proton-proton scattering but does 
not yield sufficiently precise results for proton-nucleus 
scattering at intermediate energies j23| where accuracy 
in the high partial waves is needed for convergence; in 
addition this method cannot be extended to the three- 
body problem. An improved method was also developed 
in Ref. [131 which is capable of producing more accurate 
quantitative calculations than Vincent and Phatak for 
the reaction observables for smaller values of the cutoff 
radius but still converges slowly with the cutoff radius 
for large scattering angles. An alternative method was 
proposed in Ref. |l8i] in which the limits of the cutoff 
parameter are taken analytically. More recently a novel 
technique was proposed [241 but its appHcation is still 
limited to the numerical solution of the pure Coulomb 
problem [2^ , and its numerical accuracy for high partial 
waves is yet to be tested. 

The screening and renormaHzation approach P, 01 has 
also been recently revisited leading to the treatment of 
the Coulomb interaction proposed in Ref. Q together 
with a new screening function. For completeness we 
present in here a summary of the procedure. First, we 
work with a Coulomb potential ujr^ screened around the 
separation r = R between the two charged particles. In 
this work, we choose a screening function that is different 
from the one used in Refs. J,, 3k 3\ 

i^R{r) = uj{r) [Q{R ~ r) 

+ ie(^_i?)e(3i?_^) (i + sin(|^) 



uj{r) 



(15a) 
(15b) 



where Ofe is the fine structure constant and Zp and Zt 
the ratio to the proton charge for both projectile and 
target nuclei. UnHke the screening function used be- 
fore [1, [B, 0, S, 01, this one has an analytical Fourier 
transform. In addition it possesses the same properties as 
the previous one, i.e., preserves the Coulomb interaction 
at short distances and for r > i? goes smoothly to zero. 
In Fig.[T]the shape of this new screening function is com- 
pared with the previous one, (ujnir) = uj{r)e~'''^/ ), 
for n — A. The screening radius for Eq. I|15ap corre- 
sponds approximately to the double of the former one. 
The screening radius R is chosen to be larger than the 
range of the strong interaction. However, it will be al- 
ways very small compared with the nuclear screening dis- 
tances which are of atomic scale (i.e., 10^ fm). Thus, the 
employed screened Coulomb potential ojr is unable to 
simulate the physics of nuclear screening or even model 
all features of the true Coulomb potential. However fol- 
lowing the prescription given in Ref. |l| and the techni- 
cal developments proposed in Ref. @, the results corre- 
sponding to unscreened Coulomb can be obtained. This 



procedure involves the use of a two-potential formula that 
separates the long range part from the Coulomb modified 
short range contribution. Therefore the amplitude T for 
nuclear plus Coulomb scattering reads 



lim 

R — >oo 



-1/2 



-1/2 



(16) 



where Tc is the pure Coulomb amplitude that is known 
analytically and is given in the Appendix lA II The pure 
Coulomb transition matrix that has no on-shell limit is 
not needed in the method of screening and renormaHza- 
tion P,[3,[l|. The amplitudes T^^) and T^^^ are calcu- 
lated with the nuclear plus screened Coulomb potential 
and the screened Coulomb potential alone, respectively, 
as described in Sec. II. The second term on the right side 
of Eq. (fTBl) corresponds to the Coulomb modified nuclear 
short range amplitude which is calculated numerically for 
different R, but whose R ^ oo limit is reached with high 
accuracy at finite R [S]- The renormalization factor zr 
in Eq. (fTBl) is given by 



ZR = exp (-2i</)fl), 



(17) 



where (j)R is defined in Appendix I A II 

Unlike the work in Ref. |2J, i25|] the method of screen- 
ing and renormalization has limitations at low energy 
because the screening radius has to be at least greater 
than the wavelength associated with the relative motion 
of two particles in the initial state. The lower the energy 
the higher the screening radius needed for convergence 
leading to numerical instabilities. Nevertheless for the 
typical energies of direct nuclear reactions the method 
is extremely accurate as demonstrated in the following 
section where convergence is reached at finite screening 
radius (i? < 12 fm). Although there is no a priori way 
to predict the screening radius that leads to a converged 
result, we find that convergence is readily obtained by 
monotonically increasing R until the calculated observ- 
ables change by less than a given percentage value, typi- 
cally 1% or less. 



IV. RESULTS 

In this Section we are going to apply the formaHsm 
presented above to different two-body reactions. In or- 
der to show that the method works properly for a wide 
range of energies and observables we consider first the 
calculation of p-^^Be and elastic scattering ob- 

servables obtained with a projectile-target optical model 
potential plus Coulomb at different energies {Eiat —50, 
100, 150 and 200 MeV). The results are compared with 
the solutions obtained by solving the Shrodinger equa- 
tion in configuration space, where the screening of the 
Coulomb potential is not needed. These calculations are 
performed with the code FRESCO |23| which uses the stan- 
dard method described in most Quantum Mechanics text 
books where, in each partial wave, the numerical solution 



4 



100 r r 




Figure 1: Comparison between the screening function used in 
this work with R = 5 im (full line) and the screening function 
used in Ref. Q with J? = 10 fm and n = 4 (dashed line). 



Figure 2: Differential cross section relative to Rutherford 
cross section for p-^^Be elastic scattering at 50 MeV, calcu- 
lated with different values of the screening radius R. 



of the differential equation for nuclear plus Coulomb po- 
tentials is matched, at some radius, with the appropriate 
asymptotic Coulomb wave functions which are known an- 
alytically. Convergence of results has to be tested with 
respect to the number of partial waves included and the 
matching radius . The numerical accuracy of this method 
is well known and is documented in Ref. [2^1 • 

For both reactions, we have used the parametrized op- 
tical potential of Watson j23| at the corresponding en- 
ergy per nucleon. The screening radius is taken to be 
i? = 5 fm for p-^^He, and i? = 7 fm for In order 

to perform the numerical integrations we need to intro- 
duce a certain number of mesh points as explained in 
Appendix [BJ For the p-^^Be reaction we take n^p =40, 
77.61=64, and 71^=64 and for we use n^p =50, ne=80, 

and nq=80. The CPU time needed for convergence on 
a AMD Opteron (2.4 GHz) single processor is about 10 
minutes for p-^^Be and 30 minutes for p-^^0. 

First, we show the convergence of the method itself 
with respect to the screening radius R in Fig. [51 The 
different lines are the differential elastic cross section rel- 
ative to the Rutherford cross section for the reaction p- 
^°Be at 50 MeV as the screening radius increases. From 
i? = 5 f m all curves fall on top of each other showing the 
convergence with R and the numerical stability of the 
calculation. 

Then, in Figs. [3] and [H we show the differential elastic 
cross section relative to the Rutherford cross section for 
these two reactions at the energies considered. The full 
lines represent configuration space partial-wave calcula- 
tions, and the points show the plane- wave calculations 
in momentum space. Both calculations are in very good 
agreement at all the different energies. For p-^^0 scat- 
tering we also show in Fig. [5] the analyzing power j\-y at 
the same energies as before. Again we obtain a very good 
agreement with the configuration space calculations. The 
perfect agreement between the two calculations indicates 
that the method of screening and renormalization can be 
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Figure 3: (Color online) Differential cross section relative to 
Rutherford cross section for p-^'^Be elastic scattering. Full 
lines are partial-wave calculations in configuration space and 
points are plane-wave calculations in momentum space. The 
lines a), b), c), and d) correspond to Eiat =50, 100, 150, and 
200 MeV, respectively. 



used accuratelyin high partial waves unlike the methods 
used in Refs. [li,!!!- 

Secondly we consider the reaction ^^C-^^Be. We study 
this reaction at 49.3 MeV per nucleon taking the optical 
potential used in Ref. • The screening radius is taken 
to be i? = 12 fm and the number of mesh points used are 
=80, nq=300, and ne=240 divided in three regions in 
order to have more points where they are necessary. The 
CPU time needed for convergence on the same machine 
as above is now of about 1 hour using sixteen processors. 
The correspondent elastic scattering cross section rela- 
tive to the Rutherford cross-section is shown in Fig. El 
The full line shows the configuration space partial-wave 
calculation, and the points represent the plane-wave cal- 
culation in momentum space. Both results are again in 
very good agreement. However the convergence with the 
screening radius is slower than in the previous reactions 
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Figure 4: (Color online) Differential cross section relative to 
Rutherford cross section for p-^^O elastic scattering. Full lines 
are partial-wave calculations in configuration space and points 
are plane- wave calculations in momentum space. The lines a), 
b), c), and d) correspond to Eiab =50, 100, 150, and 200 MeV, 
respectively. 
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Figure 5: Analyzing power for elastic scattering. Full 

lines are partial-wave calculations in configuration space and 
points are plane-wave calculations in momentum space. The 
lines a), b), c), and d) correspond to Eiab —50, 100, 150, and 
200 MeV, respectively. 



and slower than using a partial-wave dependent renor- 
malization factor as in Ref. 



V. SUMMARY AND CONCLUSIONS 

The two-body scattering has been studied by solving 
the Lippmann-Schwinger equation in momentum space 
without partial-wave decomposition. The Coulomb and 
spin-orbit interactions have been included. The screening 
and renormalization procedure has been used for includ- 
ing the Coulomb interaction. 

The method has been appHed to different reactions, 
p-i^Be, at 50, 100, 150, 200 MeV, and ^'^C-^°Be 

at 49.3MeV/u. The results are in good agreement with 




Figure 6: Differential cross section relative to Rutherford 
cross section for ^^C-^'^Be@49.3MeV/u elastic scattering. Full 
line is partial-wave calculation in configuration space and 
points are plane-wave calculation in momentum space. 

those obtained through the solution of configuration 
space equations using the partial-wave decomposition. 

This work shows that this procedure is reliable for a 
two-body reaction and encourages its extension to three- 
body reactions where the partial wave decomposition be- 
comes unstable at high energies or when two of the three 
particles are massive. 
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Appendix A: FOURIER TRANSFORM 

In this Appendix we present the analytical expressions 
for the Fourier transform of the potentials used in this 
paper. 

For a central potential V{r), the Fourier transform is 
a scalar function 

V{q\ q, V) = ^ j d?rV{r)e-^^^-^y\ (Al) 

where y is the cosine of the angle between q and q'. 
Integrating in the polar angles (0, (p), one gets 

which involves a single r integral. 

1. Screened Coulomb potential 

If we define the screened Coulomb potential as in 
Eq. (|15ap . the Fourier transform can be performed an- 
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alytically. Starting from Eq. (|A2p . the integral is 



UeZpZt 



1 - 



TT^ COS (2fci?) COS (fci?) 



where k = |q — q'|. When A: — > 



7r2 - (2fci?)2 



4^2 (5 - 8/7r )R . 



The renormalization phase (f)R is given in Ref. [j| 



0fl = «(9o) / ujR{r)dr, 

Jl/2qo 

which for the ti^_R,(r) given in Eq. I|15ap becomes 



(f>R = k(9o) {ln(2goi?) 



(A3) 
(A4) 

(A5) 
(A6) 



with 



and 



ln(3) + 5"^ 



Si 



(I 



S'i ( Y 1 - S'i (I) « 0.23761058, 



K{qo) = aeZpZtn/qo. 



The pure Coulomb amplitude Tc needed in Eq. I|T6|) is 
given as a function of the cm. scattering angle Oc.m. by 



(A7) 



(AS) 



^^(90) 



T{l + iK{qo)) 



2(27r)2^go sin^ (i^cm.) T (1 - ^K((7o)) 
X exp (-2iK((7o) In (sin (i6'c.„i.)))- (A9) 

2. Short-range Coulomb potential 

The Coulomb potential inside the nucleus is usually 
taken as the Coulomb potential for a uniformly charged 
sphere of radius rg whose difference from point Coulomb 
is 



2r„ I, rl 
Again starting form Eq. l|A2p we get 



^cR{q',q, y) 



aeZpZt 
" 27r2fc2 



(AlO) 



(All) 



(kro)^ + 3kro cos (fcro) — 3 sin (fcro) 
Ik^^ ■ 



When k^O 



i^cR{q',q, y) 



O-eZpZt 

207r2 



(A12) 



3. Woods-Saxon potential 

The Woods-Saxon potential, usually used for optical 
potentials, has the form 



Wo 



('-'-0) 



leading to the Fourier transform 



Vws{q',q,y) = 



Vo 



r sin (kr) 
^TT K ./g 2 + e — s — 



dr 



Defining x = r/a^ b — rg/a, and c = ka we have 



Vwsiq',q,y) 



27r2/c 



2 /"CO • / \ 

' X sm [cx) 

x — b ' 



dx- 



1 + e 



(A13) 



(A14) 



(A15) 



which, except for the constant vqo? / 2it'^ k , may be calcu- 
lated as the imaginary part of the integral , along the 
X axis. 



dx- 



1 + 



(A16) 



In order to calculate this integral we consider the integral 
in the complex plane /(z) over the first quadrant 



(2) 



dz- 



1 + 



(A17) 



that equals the sum of the same integral over the positive 
axes X and y because the integration over the arc is zero. 
Applying the Cauchy's theorem we have 



ze- 



(A18) 



n=0 



) z=b+iTT(2n+l) 



= 27re 



icb 



1 - e-2c7r 

Now we need to calculate I(^y) 



l + e-2c^ 
TT ib 

1 - e-2c^ 







dz- 



ze 



r,Z—b 



dy : — r. 

1 + e'y-'' 



(A19) 



1 + 

Using the Taylor series for the function f{y) = 1/(1 + y) 
00 ^00 



n=0 

^(-l)"e-^"/ dyye 



n^{iy-b)n 



(A20) 



-y(c~in) 



n=0 

00 



n=0 



— + 2mc 
(n2 + c2)2 ■ 
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This series converges very fast because it has a negative 
exponential increasing with n. Therefore our integral is 



for which we use the Taylor series expansion of the func- 
tion f{y) = 1/(1 + y)2 to obtain 



Im[/(^)] = Im[/(^) - /(y)] 



= 2tt 



1 - e-^c'^ 

OO 

2c^(-ire 



(A21) 

1 + e~^°'" 
7^ sin (cb) r27:r - b cos (c6) 



1-e 



-bn 



When fc ^ 



Im[/(a;) 



(A22) 



OC - 

2a5:(-l)"e-''"-3 



4. Derivative of Woods-Saxon potential 

The derivative of the Woods-Saxon potential is usually 
defined as 



dr 



(l+e(r-r-,)/a)2- 



(A23) 



Now we follow the same procedure as in the former sub- 
section but with a new function. Defining again x — r/a, 
h = Ts/a, and c — ka, we have 



r, , , , 4vsa'^ f°° , e^~^a; sin (ex) 



27r2A: 



(A24) 



from which we can define 



dz 



as the integral over the first quadrant leading to 



27riRes I — 



^f^tcz^z-b 



(A26) 



27ri ^ ({icz + l)e 



icz ^b— z^ 



z=fc+i7r(2ri+l) 



27rie' n no 

(1 - e-2c7r)2 
i6c(l - 6-2^^^)] . 



e-2='"(l + C7r) + (c7r-l) 



Now we need to calculate // 



(v) 



7^,-1 = / dz— —rz^ = / dy 



7_--(l + e-b)2 -«(l + g.y-6)2' 



(A27) 



= V(-l)"+^ne-''" / ye 



(iy-b)n 



(A28) 



-y(c-in) 



^(-l)"+ie- 



■6n 



/O 

„2 „2 



e - n + 2me 
(n2 + c2)2 • 



This series converges again very fast because it has a 
negative exponential increasing with n. Therefore our 
integral is given by 



Im[/(x)] = Im[/(z) - /(y)] 



(A29) 



= 27r- 



+ cos {cb) (e"2^''(l + e7r) + (evr - 1))] 

OO 2 

- 2e5](-l)"+ie- 



-hn 



{v? + e2)2 



When A: ^ 

Im[^(x)] 



(A30) 



OO ^ 

2a5:(-l)"+ie-^"-. 



5. Spin-orbit potential 



The spin-orbit potential, given by Eq. lfTOl) . is not cen- 
tral so Eq. (|A2p is not valid. In this case, the Fourier 
transform of the spin-orbit term is 



hx)+I{y), (A25) (q'|w,„|q) 



1 



(27r)3 



dVdVe-*i' (r)5(r' - r)e"'-' 



(A31) 



Taking into account the form of Vso (r) given in Eq. JTO] 
and the definition ofl = rxp = rx (— zV), we have 

(q>.o|q) = ^ j d^d^'e-'^ '^'vUr) (A32) 

X ,5(r' - r) [r X (-iV)] • ae^i'"". 

Integrating in r' and using the property (a x b) ■ c = 
a • (b X c) we get 

^q>-|q) = T^T / d^'r vUr)v • [q x a] e'(^-^^-\ 

(A33) 

(A34) 



(27r)3 ^ 
Since Vso{r) has the form 

Vso{r) = 



1 ljfu.a(r) 
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one may write 
leading to 

(q>so|q; 



Vsoir)r = -Vvn,s{r), 



(A35) 



1 



(27r)3 

Integrating by parts 



(A36) 



{^'\vso\q) 
-1 



j d\V ■ (v^s{r) [qx a] e'(i-i') 



(A37) 



(27r)3 

- z(q-q')-[qX(T]ydVz;„,(r)e'(i-i')-- , 
and applying Gauss' theorem, 

j (frV (v^s{r) [q x a] e'(i-<i') "-) (ASS) 
Vws{r) [q X cr] ■ dS* = 0, 



together with 

(q - q') ■ [q X a] = -q' • [q X a] - - [q' x q] • d, (A39) 
we finally have 

(q>.o|q) = ^ [q' X q] • ^1 d\ v^,{r)e^'^^-^^-\ 

(A40) 



where the integral is equals to the Woods-Saxon Fourier 
transform developed in Subsection I A 31 



Appendix B: SOLUTION METHOD 

The integral equations for the transition amplitude 
with central interaction Eq. ^ and including spin-orbit 
interaction Eq. lfT3|) are solved using the well known 
method of Fade summation j2§| with a choice of an ap- 
propriate mesh for each variable in Eq. ^ or Eq. (fT3|) . 
For the momenta q we take a Gauss-Chebyshev mesh 
converting the interval [—1,1] into [0, oo) via 



1 



(Bl) 



where b is a scale used to extend or compress the mesh. 
The scale used in this work is 5 fm~^. For the cosines x 
we take the Legendre mesh in the interval [—1, 1]. Some- 
times, depending on the problem, it is more convenient 
to divide the x-interval in regions and define a Legendre 
mesh in each one with different number of points. This 
procedure allows to increase the number of mesh points 
where they are more necessary. 

The integral in Lp in the Eq. ([7]) and lfT4|) is also calcu- 
lated with a Legendre mesh in the interval [0, 27r]. In the 
case of a central potential the integration can be done in 
the interval [0, tt] and multiplying by two. This is very 
useful if the particles do not have spin. 
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